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ABSTRACT 

We present a study of three Lya emitting galaxies (LAEs), selected via a narrow-band survey in the GOODS 
northern field, and spectroscopically confirmed to have redshifts of z ~ 5.65. Using HST ACS and Spitzer 
IRAC data, we constrain the rest-frame UV-to-optical spectral energy distributions (SEDs) of the galaxies. 
Fitting stellar population synthesis models to the observed SEDs, we find best-fit stellar populations with 
masses between ~ 10 9 - 10 10 M Q and ages between ~ 5- 100 Myr, assuming a simple starburst star formation 
history. However, stellar populations as old as 700 Myr are admissible if a constant star formation rate model is 
considered. Very deep near-IR observations may help to narrow the range of allowed models by providing extra 
constraints on the rest-frame UV spectral slope. Our narrow-band selected objects and other IRAC-detected 
Z ~ 6 /'-dropout galaxies have similar 3.6 /im magnitudes and z' — [3.6] colors, suggesting that they posses 
stellar populations of similar masses and ages. This similarity may be the result of a selection bias, since the 
IRAC-detected LAEs and /'-dropouts probably only sample the bright end of the luminosity function. On the 
other hand, our LAEs have blue /' - z' colors compared to the /'-dropouts, and would have been missed by 
the /'-dropout selection criterion. A better understanding of the overlap between the LAE and the /'-dropout 
populations is necessary in order to constrain the properties of the overall high-redshift galaxy population, such 
as the total stellar mass density at z ~ 6. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

The star formation history of high-redshift galaxies is 
an important problem in observational cosmolo gy. It pro- 
vides tests for galaxy formation the ories (e.g. lSomerville etafl 
1200 11 ISpringel & Hernauistll2003l). constrains the sources of 
reionization (e.g. Barkana & Loeb 2000, 2001), and may 
even shed light on th e nature of the "first stars" (e.g. 
iBromm & Larsonl2004l) . 

Recently, the study of star formation at high redshift has 
been catalyzed by the Infrared Ar ray Camera (IRA C) on 
board the Spitzer Space Telescope jFazio et all 12004). The 
extraordinary sensitivity of this instrument has enabled the 
detection of galaxies out to z ~ 7. At such high redshifts, 
optical and near-IR observations accessible to ground based 
telescopes and the HST sample the rest-frame UV, which pro- 
vides important information on the on-going star formation 
and dust extinction. On the other hand, IRAC samples the 
rest-frame optical light, which is less susceptible to dust ex- 
tinction and is more representative of the emission from the 
less massive stars making up the bulk of the stellar popula- 
tion. Infrared observations are therefore indispensable when 
studying the star formation history of high-redshift galaxies. 

Already, a sizable number of IRAC-detected z ~ 6 
galaxi es have been disc o vered a nd studied i n de- 
tail jDow-Hygelund etail l2005t ICharvetalJ | 2005j 
Egami et alJ 120051 lEvles et alJ 120051 iMobasher et alJ 20051 
Schaerer & Pell6ll2005t lYan et alJl2005t iMcLure et alll200ft 
Yan et al 12 006). The infrared observations provided by IRAC 
has allowed the galaxies' rest-frame UV-to-optical spectral 
energy distributions (SEDs) to be constrained. Estimates of 
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the masses and ages of the galaxies' stellar populations can 
then be obtained. A stunning result from these studies is that 
many of these high-redshift galaxies are massive (> 10 10 
M Q ), and possess stellar populations that are older than 
several hundred million years at a time when the age of the 
universe is less than 1 Gyr. 

The vast majority of the IRAC-detected z ~ 6 galaxies stud- 
ied so far have been found using the Lyman break technique 
JSteidel etalJll995l 119961 120031) extended to high redshift 
(the /' -dropout technique; Bouwens et alJ2003llStanwav et alJ 
120031 1 Yan & Windhorstll2004l) . Another, perhaps more effi- 
cient, method to search for high-redshift galaxies is to select 
candidates based on thei r Lya em ission ( Rhoads & Malhotra 
| 2001f [Aiiki et alJl2003t iHu et al]l2o5l Malhotra & Rhoadsl 
2004; Taniguchi et al. 2005). The two methods suffer from 
different selection biases. For instance, while Lya searches 
may allow discovery o f sources not detected in the contin- 
uum ( Fvnb o et al.l2001l) . they may prefer entially select young^ 
galax ies in a dust-fre e environment ( Malhotra & Rhoads 
2002). In a recent study, Gawiser et al. ( 2006) found that Lya 
emitting galaxies at z ~ 3 may be younger and less massive 
than Lyman break galaxies at similar redshifts. It is therefore 
interesting to compare the z ~ 6 Lya selected galaxies to the 
/'-dropout sample. 

In this paper, we present a study of three z ~ 5.7 galaxies 
discovere d in the Great O bservatories Origins Deep Survey 
(GOODS; Dickinson et all2003l) northern field via a narrow- 
band Lya survey. Using HST ACS and Spitzer IRAC data 
from GOODS, we study the rest-frame UV and optical prop- 
erties of Lya emitting galaxies at z ~ 5.7. The goal is to derive 
stellar mass and age estimates, and to identify differences, if 
any, compared to /'-dropout galaxies at similar redshifts. 
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FIG. 1. — Images of the three IRAC-detected z ~ 5.65 objects in i', z', 3.6, and 4.5^tm. Each panel is 10" across, except for the zoomed z' images, which are 
2" across. The images are rotated from the original mosaics, when necessary, so that north is up. The ACS stamps are binned with a 3 X 3 pixels box to bring 
out the faint extended sources. The full resolution unbinned stamps are shown in the zoomed panels. The circles indicate the apertures used in the photometry 
measurements. 



The paper is organized as follows. In § |2] we describe 
the candidate selection strategy and the photometry measure- 
ments. In § 13 we present results from population synthesis 
modeling of the observed SEDs. We compare our sample of 
Lya emitting galaxies to the /'-dropout sample in § 0] and 
the conclusions are presented in § |5] Hereafter, we will re- 
fer to the Lya emitting galaxies as Lya emitters (LAEs), and 
the /'-dropout galaxies simply as /'-dropouts. All magnitudes 
quoted are in the AB magnitude system. We adopt a cos- 
mology of {il m ,flA,h} = { 0.3,0.7,0.7|, consis tent with the 
recent results from WMAP ( Spers el et al.l2006l) . 

2. CANDIDATE SELECTION AND DATA REDUCTION 
2.1. Candidate Selection and Follow-up Spectroscopy 

The high-redshift candidates are selected using a wide-field 
narrow-band survey. At the time of this paper's preparation, 
the survey is being expanded to allow for deeper and more 
uniform selection criteria over a wider region. Details of the 
completed survey, selection criteria, and catalog will be pre- 
sented in an upcoming paper (Hu et al., in prep.). Here we will 
briefly summarize the portion of the survey and the selection 
criteria used to construct the current sample. 

The narrow-band survey, carried out using SuprimeCam 
on the Subaru Telescope, covers a 700 arcmin 2 area encom- 
passing the GOODS northern field. This data set consists of 
broad-band continuum UVBRIz' images (Cao ak etalJ l2004: 
Hu et al., in prep.), supplemented by narrow-band observa- 
tions taken with the NB816 filter (Hu et al., in prep.). The 
NB816 filter has a width of 120 A FWHM and is centered 
around 8150 A, corresponding to Lya at z ~ 5.7. 

The narrow-band observations allow for the selection of 
high-redshift candidates based on the objects' narrow-band to 
broad-band flux excess. A detailed discussion of the general 



selecti on strategy for z ~ 5.7 LAEs can be found in H u et alJ 
(2004), while considerations specific to this data set are dis- 
cussed in Hu et al. (in prep.). To summarize, objects with 
I — N > 0.7 are selected down to a narrow-band magnitude 
of N = 25.5. Due to Lya forest absorption and the galaxies' 
intrinsic continuum break, the high-redshift candidates must 
also satisfy either a) R-z' > 1.8 and be undetected in B and 
V, or b) be undetected in all passbands redward of /. 

Candidates satisfying the selection criteria are followed up 
using the DEIMOS spectrograph on the Keck Telescope. At 
the DEIMOS resolution, a number of candidates show sin- 
gle asymmetric broad lines, with long red tails and steep blue 
drop-offs. Such asymmetric lines are characteristic of Lya 
emission at high-redshift, and objects displaying such emis- 
sion lines are classified as z ~ 5.7 LAEs. 

Our initial sample of spectroscopically confirmed z ~ 5.7 
LAEs consists of 20 objects (there are additional LAEs from 
the expanded survey region), with 12 of them falling inside 
the 10' x 16'5 region where deep IRAC observations are avail- 
able from GOODS. Three objects, LAE#07, LAE#08, and 
LAE#34, show significant emission in the IRAC 3.6 and 4.5 
/im channels. From now on, our discussion will be focused 
on these three IRAC-detected z ~ 5.7 LAEs. Spectroscopic 
redshifts of the objects are tabulated in Tabled Coordinates, 
spectra, and Lya properties for the objects will be presented 
in an upcoming paper (Hu et al., in prep.) upon completion of 
the survey. 

2.2. X-Ray Nondetection 

In order to check for the possibility that our objects are 
high-redshift A GNs, we turn to the 2 M s Chandra Deep Field 
North (CDF-N; lAlexander et al.ll2003l) . in which all three of 
our objects lie. We cross-checked our objects with the point- 
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FIG. 2. — Object images in the IRAC 3.6 and 4.5 fim bands with the 
neighboring objects subtracted. 

source catalog by Al exander et al.l (T2003 ) and found that none 
of our objects is detected in that survey. The three IRAC- 
detected objects lie near the edge of CDF-N, where the sen- 
sitivity is lower at ~ 3 x 10~ 16 erg cm" 2 s" 1 (3a , full-band 
0.5 - 8 keV; see Fig. 18 in Alexan der etafl feCXB). Nonethe- 
less, a powerful AGN/quasar would have a rest-frame 0.5 - 
8 keV luminosity of arou nd 10 44 erg s" 1 (Bargeret alJ l2003l 
iBrandt & Hasingerll2005l) . which translates into an observed 
0.5 - 8 keV flux of 4 x 10" 16 erg cm -2 s _1 at z = 5.65 (as- 
suming a photon index T = 1.8). Therefore, if our sources 
were x-ray luminous quasars we would expect to detect them. 
However, we cannot rule out the possi bility of a weak AGN 
contribution. On the other hand, Wang et al. (2004) surveyed 
using Chandra a sample of 100 LAEs at z ~ 4.5 found by the 
Large Area Lyman Alpha (LALA) survey, and found no x-ray 
emission from any of the LAEs. Stacking analysis also reveals 
no detectable x-ray emission, placing a 3a upper limit on the 
average 2-8 keV luminosity at L < 2.8 x 10 42 erg s . Our 
objects are selected in a similar manner as the LALA LAEs, 
and so should be higher redshift counterparts of the LALA 
objects. Theref ore, given the objects' nondetection in CDF-N 
and the Wang et al. (2004) results, we consider it unlikely that 
our objects are high-redshift AGNs. 

2.3. Photometry 

The optical photometry is measured from the vl.O release 
of the HST ACS data from GOODS 1 A detailed descriptio n 
of the GOODS ACS data is given in Giavali sco~etall J2004I) . 
To summarize, the northern field of GOODS covers an area 
of - 160 arcmin 2 in four passbands: F435W (B), F606W (V), 
F775W (/') and F850LP (z'). The images have PSF FWHM 
~ 0."12 and reach depths of 27.4, 27.5, 26.8, and 26.5 mag 
(3cr in 1 "-diameter apertures) in the B, V, i', and z'bands, re- 



spectively. 

Infrared coverage is provided by deep Spitzer IRAC obser- 
vations, covering the same field as the ACS imaging, obtained 
as part of the GOODS program (Dickinson et al., in prep.). 2 
The data consist of imaging in four passbands, centered at 
3.6, 4.5, 5.8 and 8.0 ^m. The IRAC photometry are mea- 
sured using images from our own independent reduction of 
the GOODS data. The 3er magnitude limits (in 2"4-diameter 
apertures) in our images are 26.2, 26.4, 25.0, and 25.1 mag 
in the 3.6, 4.5, 5.8, and 8.0 /mi bands, respectively. The PSF 
FWHM is ~ 2" on average, but is slightly better at ~ 1 "8 in 
the 3.6 and 4.5 /im bands. 

The three objects are detected only in the i', zf, 3.6 /im, and 
4.5 /im bands. Stamp images of the objects in these four pass- 
bands are shown in Figure^ LAE#07 and LAE#08 appear to 
be marginally resolved in the ACS images, with approximate 
angular sizes of ~ 0."7 and 0"35, respectively. Taken at face 
value, the angular sizes correspond to physical sizes of ~ 4 
and 2 kpc for L AE#07 and LAE#08 . 

We perform aperture photometry on the ACS images for 
each of the three objects. Because of the extended nature 
of the objects in the ACS images, extra care has to be taken 
when selecting the appropriate aperture size. We adopt a 1"- 
diameter aperture in the measurements. Curve of growth anal- 
ysis of the three objects suggests that < 1% of the total flux 
fall outside an 1 "-diameter aperture. While a large aperture 
is not ideal for minimizing background noise contributions, it 
is necessary in order to measure the total flux of the objects. 
Since the fraction of enclosed flux within our chosen aper- 
ture is close to unity, and the objects are intrinsically faint and 
extended, we opted not to apply aperture corrections to the 
photometry. 

Photometry measurements on the IRAC images are consid- 
erably more challenging. At the IRAC resolution, all three ob- 
jects are unresolved and can be treated as point sources. How- 
ever, because of the large PSF, neighboring objects may sig- 
nificantly contaminate the photometry. For example, LAE#07 
is heavily blended with two other faint neighbors, while 
LAE#08 is located in the wings of a bright foreground galaxy 
(Fig.EJ. 

To address this problem, we use a deblending technique 
whereby the contaminating neighbors are subtracted using the 
high resolution ACS images convolved with the IRAC PSF. 
We use the ACS z'-band image as input since it is the closest 
to the IRAC bands in wavelength and the objects are brightest 
in this passband. In fact, we in general use the z'-band image 
as the reference for object positions. The IRAC PSF is ob- 
tained directly from the images by stacking ~ 10 bright point 
sources within the field. 

For each contaminating object in the vicinity of the target, 
the z'-band image is convolved with the IRAC PSF. Then, af- 
ter first subtracting all the neighboring sources using the cur- 
rent best-fits, the convolved image of the contaminating ob- 
ject is fitted to its IRAC counterpart. The parameters fitted 
for are the object's amplitude, background, and position. This 
process is iterated until a converged solution for every con- 
taminating object is found. Using the best-fit solution, the 
contaminating objects are subtracted from the original IRAC 
image, leaving only the target. The result of this deblending 
procedure is shown in Figure|2] It is important to note that this 

1 http://archive.stsci.edu/pub/hlsp/goods/vl/ h_goods_vl .0_rdm.html 

2 http://data.spitzer.caltech.edu/popular/goods/Documents/ 
goods_dataproducts.html 
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TABLE 1 
Measured Flux Densities 





Redshift 


i' 


i 


3.6 fim 


4.5 fim 


LAE#07 
LAE#08 
LAE#34 


5.635 
5.640 
5.671 


0.071 ±0.047 
0.070±0.045 
0.073±0.043 


0.135±0.030 
0.124±0.033 
0.091±0.025 


0.53±0.20 
0.51 ±0.20 
0.31 ±0.20 


0.41±0.17 
0.61±0.17 
0.28±0.17 



NOTE. — Flux densities are given in units of /iJy. Quoted errors are the la uncertainties. 
In the IRAC bands, the errors include systematic uncertainties coming from confusion with 
neighboring sources. The objects are undetected in B, V, R, 5.8 fim, and 8.0 fim, with 3a 
upper limits of 0.042, 0.036, 0.063, 1.0, and 1.0 fijy, respectively. 



deblending procedure works best when the objects to be sub- 
tracted are unresolved, such as the neighbors of LAE#07. The 
procedure is less successful in subtracting the bright extended 
foreground galaxy neighboring LAE#08, most likely due to 
uncertainties and spatial variations in the IRAC PSF, and pos- 
sible wavelength dependence in the foreground galaxy's mor- 
phology. Nonetheless, Figure [2] shows that the deblending 
procedure works quite well in removing contamination from 
neighboring sources. 

After the contaminating neighbors are subtracted from the 
IRAC images, we measure the objects' photometry in 2"4- 
diameter apertures, centered at the positions derived from the 
z'-band images. The aperture size is chosen to be as small 
as possible to minimize background noise and residual flux 
from nearby objects. Aperture corrections factors of ~ 0.5 
are applied to the measured fluxes. 

Uncertainties in the ACS photometry are estimated using a 
Monte-Carlo procedure. Artificial point sources, with fluxes 
equal to the objects, are inserted into the images. The posi- 
tions of the artificial sources are assigned randomly apart from 
the requirement that they are at least 1 aperture diameter away 
from detected sources. We then apply the same measurement 
procedure we used on the objects to the artificial sources. The 
resulting dispersion in the measured fluxes will serve as an 
estimate of the uncertainties in the objects' photometry. This 
procedure provides an estimate of the errors introduced by sky 
fluctuations and sky subtraction, and it takes into account the 
correlations between adjacent pixels introduced by the drizzle 
procedure used to produce the ACS mosaics. The errors de- 
rived in this way also include the confusion noise due to faint 
undetected sources. 

For the IRAC images, there is an additional source of error 
coming from the neighbor subtraction. Ideally, we would ap- 
ply a similar Monte-Carlo procedure, repeating the neighbor 
subtraction process many times on artificial sources to obtain 
the error estimates. However, the convolution-based neighbor 
subtraction process is time-consuming, and so a full Monte- 
Carlo approach is impractical. 

We therefore try to mimic the effects of nei ghbor subtrac- 
tion u sing the PSF-fitting package StarFinder (Diolaiti et al. 
2000). The errors in the source subtraction are dominated by 
image quality issues, such as variations of the PSF across the 
field, and by the properties of the individual sources (e.g. de- 
viations from being perfect point sources). The differences in 
the source subtraction using StarFinder and our convolution- 
based procedure should be small compared to the aforemen- 
tioned uncertainties. Hence, we believe that performing the 
source subtraction using StarFinder will give a reasonable 
estimate of the errors introduced by our convolution-based 
neighbor subtraction process. 

Using StarFinder, sources are detected and subtracted from 



the IRAC images. Artificial sources are inserted randomly 
into the source subtracted image and their fluxes are then mea- 
sured. Any errors arising from the imperfect source subtrac- 
tions, as well as the statistical errors coming from the back- 
ground and its subtraction, would be included in the resulting 
measured flux dispersion. The IRAC photometry errors thus 
derived are conservative since we expect our convolution- 
based subtraction to perform better than StarFinder in the case 
where the objects are not perfect point sources. Also, the arti- 
ficial sources are placed purely randomly throughout the im- 
age, so they may overlap subtracted sources. This tends to 
overestimate the uncertainties since the errors in the subtrac- 
tion are larger near the center of the sources where the flux is 
higher. The results of the flux measurements and error esti- 
mates described above are summarized in Tabled 

Finally, we supplement the ACS and IRAC imaging using 
data from the Hawaii-HDF-N project cCapa k et all2004l) . We 
attempt to measure the objects' photometry in the R -band im- 
ages from Subaru and the HK' images from the UH 2.2m tele- 
scope. However, these images are shallower than the ACS 
and IRAC data, and none of the objects are detected (in 3"- 
diameter apertures) in either R or HK' . We include the /?-band 
upper limit (27.2 mag, 3a) in our analysis for completeness. 
The HK' -band upper limit (22.7 mag, 3cr) is too high to make 
a difference in the following analysis, and is therefore omit- 
ted. 

3. SPECTRAL ENERGY DISTRIBUTION AND STELLAR 
POPULATION 

3.1. Population Synthesis Modeling 

After we have obtained reliable photometry for the galax- 
ies, we use the stellar population synthesis model of 
iBruzual & Chariot! (120031 hereafter BC03) to model the 
galaxies' SEDs. The basic strategy is to calculate model 
SEDs, varying the mass, age, dust reddening, metallicity, and 
star formation history. The model SEDs are redshifted and in- 
tegrated through the filter response function. The model pre- 
dictions are then compared to the observed photometry and 
the best-fitting model is found by \ 2 minimization. 

For simplicity, and to facilitate comparison with other 
wo rks in the literat ure, we will assume throughout this study 
the Salpeter ( 1955) initial mass function (IMF), implemented 
in the BC03 models with lower and upper mass cutoffs of 0. 1 
and 100 M Q . We will discuss the effects of using a different 
IMF and other variations in the population synthesis models 
in § 13.31 We explore models with five different metallicities: 
Z/Z = 0.005, 0.02, 0.2, 0.4, and 1. Note that metallicity is 
not treated as a free parameter in the fitting process, but rather 
an input assumption of the model. The effect of dust red- 
dening is included using the Calzetti et al. (2000) model, with 
E(B-V) varying between 0.0 - 1.0. The effect of Lya for- 
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FIG. 3. — Best fit SSP models to the three IRAC-detected z ~ 5.7 LAEs. 
The blue-dotted (red-dashed) line is for a model with Z = Zq (0.005 Zq). 
The observed SEDs are plotted as diamonds with error bars. The horizontal 
error bars indicate approximately the FWHM of the passbands. When an 
object is undetected in a certain passband, the corresponding 3<r upper limit 
is shown. Note that a 30% fractional error is added to the /'-band to account 
for the Lya line contribution. 

est attenuation is included using the model of M adaul i 19951) . 
There is no need to fit for the objects' redshifts, since these 
are known from spectroscopy. An additional constraint aris- 
ing from the known redshifts is that the inferred ages should 
not exceed the age of the universe, which is approximately 
0.99 Gyr at z ~ 5.65. We therefore impose an upper limit of 

0. 9 Gyr on the age of the best-fit models. 

The star formation history (SFH) is parameterized by the 
time evolution of the star formation rate (SFR). We consider 
two SFHs in our analysis: Simple Stellar Population (SSP), 

1. e. instantaneous burst, and constant SFR. We choose not to 
consider more complicated SFHs, such as an exponentially 
declining SFR, because our ability to constrain extra parame- 
ters is limited. All our objects are detected in only four pass- 
bands, so we can constrain up to three parameters, which we 
choose to be mass, age, and dust reddening E(B-V). Using 
an exponential SFR would introduce an extra parameter, the 
star formation timescale t, which can be constrained only if 
we fix the value of another par amet er such as E(B-V). We 
will return to this point later in § 13.31 Note that the SFHs used 
in our analysis (SSP and constant SFR) represent the two ex- 
tremes in star formation timescales, and an exponential SFR 
can be thought of as an intermediate between these two exam- 
ples. 

With two SFHs and five metallicities, there are 10 differ- 
ent classes of models in total. The best-fitting model within 
each class is found by x 2 minimization. There are three free 
parameters (mass, age, and reddening) and four data points 
(if, z', 3.6 /im, and 4.5 fim), resulting in 1 degree of freedom 
in the x 2 fits. In practice, a commonly used measure of the 



FIG. 4. — Best fit constant SFR models to the three IRAC-detected z ~ 5.7 
LAEs. Symbols and line styles have the same meaning as in Figure l3l 

goodness of fit is xt< the x 2 per degree of freedom. When 
v = 1, a model is rejected at the 1, 2, and 3-cr levels if xt is 
larger than 1, 4, and 9, respectively. We in general consider a 
fit to be acceptable if it gives xt < 4- In certain cases when 
v = 2 (see § I3.2> . a fit is deemed acceptable if xl < 3.09 (the 
corresponding 2a level for v = 2). 

In order to make full use of the data available, we incor- 
porate undetected data points as follows. Given the detection 
limit of the exposure a and the model flux in the channel /i, 
we calculate the probability I = P(< er|/z) that the actual flux 
of the object is less than a. For the purpose of this calcula- 
tion, the flux is assumed to be normally distributed with mean 
[i and width a. For each undetected channel, a value of-ln(Z) 
is then added to the total x 2 - in this way, if the model flux is 
much higher than the detection limit, /i ^ a, I would be small 
and the model would incur a large x 2 penalty. For example, 
if the model flux in a certain channel is at the 3cr detection 
limit, /x = 3ct, the x 2 penalty would be 3.8. In practice, the 
undetected data points have only a small effect on the total 
X 2 , since the detection limits in most channels can easily be 
satisfied by reasonable models. 

Since the Lya emission line is located inside the /'-band 
at z ~ 6, and the BC03 model does not predict emission line 
strengths, we have to take into account the possible Lya line 
emission when fitting to the data. The fractional contribu- 
tion of the line to the broad-band flux can be approximated as 
F L /F B = (R- l)W N /W B , where W N /W B and R = F V {N)/F V {B) 
are the ratios of the narrow-band to broad-band widths and 
flux densities, respectively. Our objects are selected based 
on their narrow-band to broad-band flux excess, and typically 
have N—I ~ -1.5. Substituting W N « 120 A and W B « 1400 
A for the NB816 and / filters, and R = lO" - 4 ^" 7 ' « 4, we find 
that the Lya line may contribute around 30% of the total I- 
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FIG. 5.— A\ 2 contour for LAE#07 in the Z SSP model. The diamond 
marks the best fitting model. The solid, long-dashed, and short-dashed lines 
represent the 1, 2, and 3-ct confidence regions. 

band flux. Note that this approximation implicitly assumes 
that the Lya forest truncates similar fractions of the narrow- 
band and broad-band fluxes. Since Lya at z ~ 5.7 is located 
slightly redward of the /'-band center, Lya forest truncation 
affects the /'-band flux more severely than / or NB816. The 
Lya line contribution may therefore be substantially higher 
than 30% in the /'-band. Nevertheless, we add an optimistic 
30% fractional error to the /'-band flux when we perform the 
X 2 fits. 

3.2. Fitting Results 

The observed SEDs of the three galaxies are shown in Fig- 
ure|3] Ignoring the lines for the moment, the first thing to take 
note is that all three galaxies exhibit similar rest-frame UV 
and optical emission properties. This is not surprising since 
the three objects are selected from same data set using the 
same selection criteria. The galaxies' SEDs are characterized 
by red UV to optical colors, with the flux increasing by more 
than a factor of 3 from z' to 3.6 /im. The SEDs also show 
signs of flattening out beyond 3.6 /im. 

The best fit BC03 models to the observed SEDs are shown 
in Figure[3]and Figure|4]for the SSP and constant SFR mod- 
els, respectively. Only models with the two extreme metal- 
licities (Z Q and 0.005 Z ) are shown, even though fits were 
performed for all five possible metallicities. In practice, we 
find good fits with competitive x 2 for all five metallicities. 
Note that the /'-band flux is consistently higher than the mod- 
els, owing to the Lya line emission. Also, since the /'-band 
straddles the Lyman-break, the model band-integrated flux is 
larger than it seems, and is usually within la of the observed 
value. 

The best-fit parameters differ depending on the SFH and 
metallicity (Table |2ji. The best-fit Z Q SSP models to the 
three LAEs have ages ~ 5 Myr, masses ~ 5 x 10 9 M Q , and 
E(B-V) ~ 0.3-0.4. The best-fit 0.005 Z Q SSP models favor 
slightly older ages (~ 50 Myr), larger masses (~ 10 10 Mq), 
and somewhat lower extinction. This is mostly because metal- 
poor stars produce more UV photons, so the best-fit models 
need to be older in order to match the red UV-to-optical colors 
in the data. These results suggest that our sample of z ~ 5.7 
LAEs displays qualities of dusty young galaxies, seen imme- 
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FIG. 6.— Best fit SSP models to the three IRAC-detected z ~ 5.7 LAEs 
with fixed E(B — V) = 0. Symbols and line styles have the same meaning as 
in Figurel3l 

diately after or during a burst of star formation. 

Compared to the SSP models, the constant SFR models 
in general yield older ages at > 700 Myr and slightly larger 
masses at a few times 10 ll) M Q . In some cases, however, the 
best-fit constant SFR models have similar parameter values as 
the SSP models. This is not surprising since for young ages 
the constant SFR model is similar to the SSP model. 

Given that we only have four data points, some of which 
have fairly large error bars, it is important to explore the range 
of models allowed by the data. Furthermore, the somewhat 
high levels of extinction we obtained may seem be a little sur- 
prising. Several previous studies have fo und littl e or no dust 
extinction in z ~ 6 /'-dropouts Jp ow-Hygelund et a^ |2005t 
lEvles et alJ2005llMobasher et alJ2005tlYan et all2005l) . The 
fact that our galaxies are selected based on their Lya emission 
would also lead one to expect low extinction va lues, since Lya 
photons are very susceptible to dust scattering dCharlot & Fall 
Il99ll IChen & Neufeldll99i lHansen & Ohll2006h . Do the 
data permit models with low extinction? What is the range 
of age and mass allowed by the data? 

To answer these questions, we calculate confidence regions 
for the best-fit parameters. In Figure |5j we show the A\ 2 
contours for LAE#07 in the Z Q SSP model. It can be seen 
from the figure that firm upper limits on the age and mass 
can be placed on LAE#07. In the context of the SSP model, 
LAE#07 has a mass of < 3 x 10 10 M Q , and an age of < 300 
Myr. The contours and limits for the other two objects are 
similar. 

One important thing to note in Figure [5] is that there is a 
range of models that fits the data. This is the result of a funda- 
mental degeneracy in the models: the red UV-to-optical colors 
of the observed SEDs can be satisfied by either a significant 
Balmer break (implying a relatively mature stellar popula- 
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TABLE 2 
Best-Fit Model Parameters 





Age [Myr] 


Mass [Mq] E(B-V) 


2 

X !J 


Age [Myr] 


Mass [Mq] E(B-V) 


2 

A V 






SSP, Zq 






SSP, 0.005 Zq 




LAE#07 


4.8 


2.4x10 s 0.275 


2.05 


60 


1.7xl0 lu 0.200 


1.97 


LAE#08 


3.2 


5.0xl0 9 0.425 


2.19 


30 


2.9 XlO 10 0.350 


2.20 


LAE#34 


4.4 


1.4X10 9 0.275 


2.42 


80 


9.1x10 s 0.150 


2.36 






Constant SFR, Zq 






Constant SFR, 0.005 Zq 




LAE#07 


5.0 


4.5 x10 s 0.375 


2.14 


720 


2.6xl0 lu 0.175 


2.19 


LAE#08 


4.8 


6.9 xlO 9 0.425 


2.19 
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3.9xl0 10 0.225 


2.18 


LAE#34 


720 


1.1 xlO 10 0.100 


2.47 
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1.7xl0 10 0.150 


2.37 






SSP, Zq, No Extinction 




SSP, 0.005 Z Q , No Extinction 




LAE#07 


90 


9.0x10 s 0.000 


1.25 


260 


1.7xl0 lu 0.000 


1.05 


LAE#08 


130 


1.4xl0 10 0.000 


1.45 


320 


2.3 xlO 10 0.000 


1.28 


LAE#34 


80 


5.1x10 s 0.000 


1.41 


230 


9.2 x10 s 0.000 


1.27 



tion), or a young stellar population with high extinction. Fig- 
ure |5] shows that even though the best-fit model has a young 
age and high extinction, there are models with lower extinc- 
tion values and older ages that will fit the data equally well. 

Because of this degeneracy, we performed an alternative fit 
to the data forcing E(B - V) = 0. The results of this fit for the 
SSP model are presented in Figure |6] Good fits to the data 
can be obtained for all metallicities. The best-fit models to 
the three galaxies in the Z Q case have masses ~ 10 10 M Q , 
and ages ~ 100 Myr. The 0.005 Z Q model gives older ages 
(around a few hundred Myr) and slightly larger masses. In 
general, the no extinction models have a lower xl ( see Ta- 
ble |2j, mostly due to the extra degree of freedom {v = 2 in 
this case). Note that in the SSP models presented here, there 
will not be enough young massive stars around to produce 
significant Lya emission becuase of the older ages. The re- 
sults we obtained should therefore be regarded as a fit to the 
older stellar population that dominates the total mass of the 
system. The observed Lya emission can still be explained by 
the presence of a young but significantly less massive stellar 
population within the galaxy. We find no satisfactory constant 
SFR models with E(B-V) = 0. Without dust extinction, these 
models produce too much UV compar e d to t he data. 

In a recent study, iLe Delliou et alJ (|2006) made detailed 
predictions of the properties of LAEs, based on a hierarchi- 
cal galaxy formation model. We find that the best-fit stellar 
masses of our LAEs are ab out 1-2 orders of magnitude larger 
that the values predicted by Le Delliou et al. ( 2006). The most 
probab le reason for this discrep ancy is that a top heavy IMF is 
used in Le Delliou et a resulting in a lower M/L ratio 

and hence lower overall mass. By requiring that our LAEs be 
detected by IRAC, we may also be sele cting galaxies from 
the high-mass end of the distribution (§ I4.2l i. Other prop- 
erties of the LAEs, such as the broad-band magnitudes as a 
function Lyq flux, agre e quite well with the predictions of the 
ILe Delliou eTal] J2006I) model. 

In summary, we find that the three IRAC-detected z ~ 5.7 
LAEs in our sample show similar emission properties and can 
be fitted with similar model SEDs. The current data are unable 
to distinguish between models with different metallicities, as 
best-fit models with competitive \ 2 can be found for all metal- 
licities. The observed SEDs are broadly consistent with two 
classes of models. In one model, the galaxies are dusty, young 
(ages < 100 Myr), and relatively less massive (M ~ 10 9 Mq). 
In the other model, the galaxies are less dusty, more massive 
(M ~ 10 10 M Q ), and older (> 100 Myr) with a significant 



Balmer break. Both types of models can provide satisfactory 
fits to the data. This degeneracy stems from the fact that we 
only have two data points (/' and z') with which to measure 
the UV spectral slope, which is sensitive to both dust extinc- 
tion and metallicity. Very deep near-IR observations in the J, 
H, or K bands would provide valuable constraints on the UV 
spectral slope, and may help to break the degeneracy and rule 
out some of the possible models. 

3.3. Alternative Models 

In addition to the basic models discussed in the previous 
section, we also tried to fit a number of alternative models to 
the data, in order to explore several effects not included in 
the basic fits. These alternative models are discussed in turn 
below. 

Exponential SFR. — We did not consider an exponential 
SFR in our basic fits because with only four data points, it 
is not possible to fit for the mass, age, reddening, and star 
formation timescale (t) simultaneously. However, if we are 
willing to sacrifice one degree of freedom in another param- 
eter, we can attempt to fit an exponential SFR model to the 
data. In the previous section, we found that the galaxies can 
be fitted by models with very little dust extinction. Other 
studies of galax ies at similar redshifts also found littl e or no 
dust extinction tD ow-Hygelund et a ll2005tlEvles et alJ2005t 
Moba sher et alJl2005HYan et alJl2005h . We therefore attempt 
to fit an exponential SFR model with E(B-V) fixed at zero. 
Good fits to the data can be found using all metallicities, with 
t on the order of 400 Myr. Since the best-fit r is fairly large, 
the exponential model is not unlike the Constant SFR model. 
In fact, the two SFH models give similar masses and ages, 
which are ~ 10 10 M Q and ~ 700 Myr, respectively. 

Chabrier IM F. — In addition to the Salpeter IMF, the 
IChabrierl ll2b03l IMF is also implemented in the BC03 mod- 
els. The Chabrier IMF has a flatter distribution than Salpeter 
at < 1 M Q . Stars more massive than 1 M Q therefore have a 
larger relative contribution to the total luminosity, resulting in 
a lower M/L ratio. Hence, the Chabrier IMF in general pro- 
duces best-fit stellar masses that are ~ 30% -40% lower than 
the Salpeter IMF. 

Alternative Population Synthesis Model. — An alterna- 
tive popu l ation synthesis model was recently introduced by 
Maraston (2005). The main new ingredient of this model is 
an increased emphasis on the contributions of thermally pul- 
sating asymptotic giant branch (TP-AGB) stars to the overall 
luminosity. Contributions from TP-AGB stars are strongest 
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FIG. 7. — z' — [3.6] vs. i' — z' colors of the LAEs in this paper and other 
IRAC-detected z ~ 6 (' -dropouts in the literature. Only the eiror bars of the 
object with the faintest z' magnitude in each group are shown. This is done to 
avoid clutter and to provide a conservative estimate of the errors associated 
with each group of measurements. The vertical dotted line corresponds to 
i' — z' = 1.5, a typical selection criterion for the z ~ 6 i' -dropouts. The curve 
shows the color evolution of a dust-free, 0.005 Zq, SSP mode. The curve 
corresponds to model ages from ~ 30 Myr near the bottom of the plot to 
~ 800 Myr near the top. 

in stellar populations around ~ 1 Gyr old, and their ef- 
fect is to raise the luminosity at near-IR and longer wave- 
lengths. In general, this results in a lower M/L ratio, and 
hence smaller best-fit masses and ages ( Mobashe r et al.l2005l 
Maras toii et al.ll2006h . However, since the LAEs in our sam- 
ple are in general much younger than 1 Gyr, and we have 
data available only in the rest-fram e optical an d shorter wave- 
lengths, we do not expect the Maraston (2005) model to give 
significantly different results than the BC03 model. On the 
other hand, TP-AGB stars could potentially contribute to the 
5.8 and 8.0 flux of z ~ 6 galaxies, and we stress that it is 
important to investigate this alternative model in cases where 
5.8 and 8.0 /im data are available. 

4. COMPARISON TO THE ('-DROPOUT SAMPLE 

One of the main goals of this study is to compare and con- 
trast the properties of the z ~ 6 IRAC-detected LAEs and i'- 
dropouts. As we have mentioned in §0 the different selec- 
tion criteria for LAEs and /'-dropouts may imply fundamen- 
tal physical differences between these two populations. One 
thing we try to accomplish is to compare the two populations 
using model-independent empirical indicators of the galaxies' 
properties, such as the rest-frame UV-to-optical color. Our 
sample of IRA C-detected z ~ 6 /'-dropouts is drawn fr om the 
works oflDow-Hygelund etlfll2005t . lEvles et al.ld2005l) . and 
lYan et al.l (120051) . Other objects are not included either be- 
cause they lie at higher redshifts (z > 6.5) so accurate i'—z 1 
color cannot be measured owing to Lya forest attenuation, or 
because IRAC photometry is not available. 

4.1. i' — z' color 

The z' — [3.6] color is plotted against the i'—z! color in Fig- 
ure[7]for both our LAEs and the z ~ 6 /'-dropouts. The most 
striking feature in Figure 0is that the LAEs have much bluer 
i'—zf color than the /'-dropouts. In fact, all three LAEs in our 
sample have /' - zf colors that are well blueward of the typi- 
cal /'-dropout selection criterion of i'—z' > 1.5. One possible 
reason for this difference is that the /'-dropouts are preferen- 
tially selected to have red i'—z! color, which may also result in 
the /'-dropouts being selected at slightly higher redshifts than 
the LAEs. Another reason is t he u nusually high /'-band flux 
observed in our sample (c.f. § I3.2> . which may be explained 
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FIG. 8. — IRAC 3.6 fim magnitude vs. — £ color of the LAEs in this paper 
and other IRAC-detected z ~ 6 £' -dropouts in the literature. 

by the Lya line contribution. We have remarked in § 13. ll that 
the Lya line can contr ibute ~ 30 % or more of the total broad- 
band flux. In addition, iHuetal] 12004) found that, in a sample 
of LAEs selected in a similar manner as our objects, there is 
a clear correlation between the narrow-band and /-band mag- 
nitudes, and that all of the spectroscopically identified LAEs 
in the sample have I—z?<l .5 (see their Fig. 5). Similarly, we 
find that all three of the LAEs have i'-z' < 1.5. 

The blue /' - z' colors of the LAEs in our sample implies 
that they would have been overlooked by previous surveys 
based on the /'-dropout selection criterion. It is therefore im- 
portant to understand the overlap between the LAE and /'- 
dropout galaxies, so that the properties of the overall z ~ 6 
galaxy population, such as the total stellar mass density, can 
be constrained. 

4.2. Stellar Population 

At z ~ 6, the zf and IRAC 3.6 /im bands sample the rest- 
frame UV and optical light, respectively. If we ignore the 
effects of dust extinction, the z' magnitude is sensitive to the 
massive young stars in the galaxy, while the 3.6 /im magni- 
tude measures the contributions from less massive, evolved 
stars. The z' - [3.6] color, plotted against the i'—z! color in 
Figure^] is therefore an indicator of the age of the stellar pop- 
ulation. The model color evolution, also plotted in Figure 
help to illustrate this point. The curve starts at ~ 30 Myr near 
the bottom of the plot. As the age increases, z! — [3.6] also 
increases, and eventually reaching z' - [3.6] = 2.5 at an age of 
-800 Myr. 

As we have mentioned before, the presence of Lya emis- 
sion may signal a young stellar population. However, Fig- 
ure shows that there is no significant difference between 
the z' — [3.6] color of the LAEs and that of the /'-dropouts. 
Therefore, we find no evidence to suggest that there is a sys- 
tematic age difference between the IRAC-detected LAE and 
/'-dropout populations. However, we stress again that the 
Z 1 — [3.6] color is a good indicator of stellar population age 
only in the absence of dust. 

Figure[8]shows the 3.6 /im band magnitude plotted against 
i'-z'- All objects in Figure [8] lie at similar redshifts, and 
hence similar luminosity distances. This implies the 3.6 /tm 
magnitude may serve as a measure of the total evolved stellar 
mass in the galaxies. The 3.6 /im magnitudes of the LAEs 
fall comfortably inside the range of magnitudes observed for 
the /'-dropouts. Therefore, we again find no evidence that the 
IRAC-detected LAEs and /'-dropouts are systematically dif- 
ferent. 
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In terms of stellar population synthesis modeling, the re- 
sults are consistent as well. For the z ~ 6 /'-dropouts, stellar 
population synthesis modeling in general gives ages around a 
few hundred Myr, and masses around 10 10 M Q . This is mostly 
in line with the values we obtained for the LAEs (Tabled, at 
least in the dust-free SSP case. 

It is important to keep in mind the comparison presented 
here is restricted to the IRAC-detected LAEs and /'-dropouts. 
By requiring IRAC detection for the galaxies, the selection 
is biased towards the bright end of the luminosity function. 
Recent studies of /'-dropouts at z ~ 6 suggest that the IRAC- 
invisible galaxies are in general younger and less m assive 
JEvles et al.l2006HYan et all20 06). Gawiser et al. <200<J also 
found that LAEs at z = 3.1 typically have lower masses at 
~ 5 x 10 8 Mq. The similarities seen in FigureQand Figure[8] 
may therefore be a result of this IRAC-selection bias. 

5. CONCLUSION 

In this paper, we studied in detail the properties of three 
Lya emitting galaxies, each spectroscopically confirmed to 
lie at z ~ 5.65. Using ACS and IRAC data from GOODS, we 
measured the galaxies' SEDs in the rest-frame UV through 
optical. Stellar population synthesis modeling then allows us 
to place constraints on the galaxies' masses and ages. Our 
main conclusions may be summarized as follows. 

The three IRAC-detected LAEs in our sample exhibit simi- 
lar emission properties, characterized by a red rest-frame UV- 
to-optical color, with the flux increasing by more than a fac- 
tor of 3 from z' to 3.6 /im. This feature can be explained 
by a young stellar population with significant dust extinction, 
or an older stellar population that has developed a B aimer 
break. Assuming the SSP model for the star formation his- 
tory, we find best-fit stellar populations with masses between 
- 10 9 - 10 10 M and ages between -5-100 Myr. How- 
ever, stellar populations as old as 700 Myr are admissible if 
a constant SFR model is considered. Very deep near-IR ob- 
servations may help to narrow the range of allowed models 
by providing extra constraints on the rest-frame UV spectral 
slope. The available data provide very little constraints on the 



LAEs' metallicity. We fitted stellar population synthesis mod- 
els using five metallicities ranging from 0.005 Z Q to Z Q , and 
found that equally good fits with similar parameter values can 
be obtained for all metallicities. 

In comparison with other IRAC-detected z ~ 6 galaxies se- 
lected based on the /'-dropout technique, we find that the 
LAEs and /'-dropouts possess similar z' — [3.6] colors, sug- 
gesting that they are similar in ages. Also, the LAEs and /'- 
dropouts have comparable 3.6 /jm magnitudes, which imply 
they have similar masses. On the other hand, the compari- 
son is restricted to the IRAC-detected LAEs and /'-dropouts, 
which are likely to be the brightest and most massive mem- 
bers of their respective populations. The observed similarities 
between the IRAC-detected LAEs and /'-dropouts may be a 
result of this selection bias. 

Even though the IRAC-detected LAEs and /'-dropouts 
share some common characteristics, the LAEs have much 
bluer i' — zf colors. Many previous searches for z — 6 galax- 
ies were based on a combination of i'—zf color selection and 
non-detections in bands blueward of /'. Because of their blue 
i'—zf colors, LAEs would be overlooked by these searches un- 
less additional selection criteria are incorporated. One present 
challenge is therefore to understand the overlap between the 
LAE and /'-dropout populations. The solution to this problem 
will shed light on the z — 6 galaxy population in general, and 
help to constrain the total stellar mass density, as well as the 
total contribution of massive galaxies to the ionization back- 
ground at z — 6. 
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